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We investigate theoretically the current induced spin-orbit torques in Ill-V diluted magnetic 
semiconductors, by combining a six-band Kohn-Luttinger Hamiltonian and a linear Dresselhaus 
spin-orbit coupling. The nonequilibrium spin density and torque are calculated using a Boltzmann 
equation taking into account transport scattering times obtained from a first-order Born approx- 
imation. For two different materials and a wide range of parameters, we scrutinize the angular 
dependence of the spin torque. We find (1) the spin torque develops a nonlinear dependence on the 
exchange coupling; (2) a strong correlation between the angular dependence of the torque and the 
anisotropy of the Fermi surface. 
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I. INTRODUCTION 

In pursuing efficient mechanisms for magnetization 
switching, spin-transfer torque attracts much attention 
in the field of spintronics. The Slonczewski-Berger pro- 
posal of spin transfer torque adopts a multilayer geome- 
try that comprises two ferromagnetic layers separated by 
a normal metal spacer where the spin angular momen- 
tum is transferred by the carriers from one ferromagnet 
(polarizer) to the other (free layer), thus inducing exci- 
tation on the latter.'^ The spin transfer torque has been 
achieved in various magnetic structures such as metal- 
lic spin valves,* magnetic metal tunnel junctions,'^ and 
especially semiconductors i^ii 

An alternative approach to accomplish magnetiza- 
tion switching uses spin-orbit interaction that couples 
orbital angular momentum to spins (with s denoting 
the spin operator) of itinerant charge carriersi^ The 
spin-orbit Hamiltonian can be generally described by 
Hso = s ■ Bso{k), where the effective magnetic field 
Bso{k) depends on momentum. Two major types of 
spin-orbit coupling, Rashba and Dresselhaus, exist in 
a broad range of materials and hybrid structures: the 
Rashba-type coupling often finds its existence in multi- 
layer structures without inversion symmetry;^ the cubic 
Dresselhaus type {Bso{k) oc k^) emerges in zinc-blende 
crystals, and strain usually introduces a linear Dres- 
selhaus (Sso(fc) oc k) spin-orbit coupling. In general, 
in a conductor where both ferromagnetic exchange and 
spin-orbit couplings are present, the nonequilibrium spin 
density generated (through the spin-orbit coupling) by 
the injected charge current is not fully aligned to the 
direction of the exchange field, thus giving rise to a 
torque, the so-called spin-orbit torque, exerted on the 
magnetizationJ^^— Therefore, a separate polarizer is not 
needed within such a scheme. Several experiments on 
magnetization switching in metallic thin films have pro- 
vided strong indication on the spin torque induced by a 
Rashba-type spin-orbit coupling, see for example, Refflg. 

Besides metallic systems, diluted magnetic semicon- 
ductor (DMS) is an excellent candidate for various 



spin-torque based applications. DMS ideally integrates, 
in a single realization, itinerant charge carriers, tun- 
able magnetism and spin-orbit couplings. Recent ex- 
periments have shown that magnetization can be re- 
versed by an electric current in strained (Ga,Mn)As 
films fi^ii^ employing exactly the spin-orbit torque. The 
magnetism of a DMS depends on magnetic anisotropy, 
which is tunable by applying an electric field, doping or 
strain engineering. I2r~ Electric fields alter the magnetic 
anisotropy by changing the carrier concentration.— Dop- 
ing affects straini^ and the exchange interaction)^ Ex- 
periments in (Ga,Mn)As have shown the amplitude of 
the spin torque increases by doping phosphor. On one 
hand, doping in general alters lattice constants, resulting 
in a change in exchange interactionj^ On the other hand, 
doping modifies band structures)^ leading to a variation 
in magnetic anisotropy. Meanwhile, by engineering the 
strain, the direction of the easy axis of the anisotropy 
field can be rotated from being in-plane to out-of-planeji^ 
thus offering an opportunity to further manipulate the 
magnetism. A tunable bulk inversion asymmetry is an- 
other characteristic feature of DMS, which allows gener- 
ation of (linear) Dresselhaus spin-orbit couplings. Spin- 
tronic properties and versatility that DMS have demon- 
strated in numerous applications^^ make them an out- 
standing platform to investigate physics related to spin- 
orbit torque. 

Despite a few studies on a similar topic ^li^i a theo- 
retical understanding of the current-induced spin torque 
in a DMS is far from comprehensive, which motives the 
present work. In this article, we focus on several issues 
that have not been addressed in existing literature and 
organize the rest of the article as follows. Section |TT] out- 
lines the model and methods on which the due discus- 
sion is based. In Sec lIIIl we discuss a strong correlation 
between an anisotropic Fermi surface and the angular 
dependence of the torque. Section HVl contributes to ef- 
fect of the exchange interaction on the spin torque. We 
highlight not only on the role of the exchange coupling 
as a mediator between the local and itinerant spin an- 
gular momenta, but also on its contribution to the spin- 
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dependent scattering. In Sec|V] we investigate the spin 
torque as a function of hole concentration. In this paper, 
we provide a comparative study using two popular DMS 
candidates: (Ga,Mn)As and (In,Mn)As. Finally, SeclVT 
concludes the article. 



II. HAMILTONIANS AND METHODS 

The system under investigation is a uniformly mag- 
netized single domain DMS film made of, for example 
(Ga,Mn)As or (In,Mn)As. We assume the system is well 
below its critical temperature. An electric field is ap- 
plied along the x direction to bring the system out of 
equilibrium. It is worth pointing out: we consider here 
a large-enough system to allow us disregard any effects 
arising due to boundaries and confinement. 

We use a six-hand Kohn-Luttinger Hamiltonian to de- 
scribe the band structure of the DMS,^^ 
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Hamiltonian Eq. ([1]) comprises subspaces of heavy holes 
(hh), light holes (Ih), and spin-orbit splitted bands (so), 
which are defined as 
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The other parameters appearing in Eq.([T]) are defined as 
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where the phenomenological Luttinger parameters 71,2.3 
determine the band structures and the effective masses 
of valence-band holes. Particularly, 73 is associated with 
the anisotropy of energy band structure around the F 
point at 72 ^ 73)^ 

Setting 71 = 1.0 while 72 = 73 = 0, the six-band 
Hamiltonian reduces to a simple free-electron model. If 
72 = 73 7^ 1 the Fermi surfaces of both minority and ma- 
jority hole bands become spherical. As we will discuss in 



Sec mil the angular dependence of the spin torque de- 
velops an interesting correlation with the Fermi surface 
anisotropy. 

The bulk inversion asymmetry allows us to augment 
the Kohn-Luttinger Hamiltonian by a strain-induced 
spin-orbit coupling of the Dresselhaus typej^i^^ We as- 
sume the growth direction of (Ga,Mn)As is directed along 
the z-axis, two easy axes are pointed at the x and y, 
respectively,-® In this case, the components of the strain 
tensor e^x and eyy are identical. Consequently, we may 
have a linear Dresselhaus spin-orbit coupling 

■f^DSOC = fi{uxkx - <yyky), (8) 

given j3 the coupling constant that is a function of the 
axial strain , ^'^•^^ '^x{y) is the Pauli matrix and fca:(i/) are 
the corresponding components of the wave vector. The 
spin-orbit coupling Eq.® can be viewed as arising from 
a momentum-dependent effective magnetic field Bd that 
is pointing at the direction {—kx,ky). 

In the DMS systems discussed here, we incorporate a 
mean-field like exchange coupling to enable the spin an- 
gular momentum transfer between the hole spin (s) and 
the localized (d-electron) magnetic moment of ionized 



Mn^''" acceptors 



JpdNunSai^ ■ S (9) 

where Jpd is the coupling constant. In this article, we 
keep Jpd positive, for the exchange coupling is believed to 
be antiferromagnetic in (Ga,Mn)As^- and (In,Mn)Asi^ 
Here Sa = 5/2 is the spin of the acceptors. The hole 
spin operator, in the present six-band model, is a 6 x 6 
matrixes! The concentration of the ordered local Mn-^+ 
moments A^Mn — 4,x/a'^ is given as a function of x that 
defines the doping concentration of Mn ion. a is the 
lattice constant. 

Therefore, the entire system is described by the total 
Hamiltonian^, 



sys 



Hkl + Hex + Hy)SOC- 



(10) 



The applied electric field is not explicitly included in 
Eq. PT])) but to be treated in the framework of linear re- 
sponse. 

In order to calculate the spin torque, we proceed to 
determine the nonequilibrium spin densities in response 
to an external electric field. In the absence of exchange 
couplings, such an electrically induced spin polarization 
has been discussed by Edelstein in an infinite medium;^ 
which is often referred in literature as the inverse spin- 
galvanic effect. A comprehensive analytical derivation of 
the spin-orbit spin torque using field theory technique 
or Kubo formalism is possible^i but beyond the scope 
of this article. In the linear response regime, by solving 
the Boltzmann equation and retaining the leading order 
correction in distribution function which is linear in the 
electric field Ex, we obtain the nonequilibrium spin den- 
sity of holes S,^ 



S = eEx^T^ 
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The expectation values of the spin and velocity operators 
are 



(s) = {lpnk\s\ilnk) 
(v) = {tpnk\v\i^nk) 



(12) 



The velocity operator is obtained using the Heisenberg 
equation v — {i/h)[r,Hsya]- Diagonalization of the total 
Hamiltonian [Eg. dTUl) ] generates the eigenstates tpnk- In- 
dex n labels 6 different bands: spin-splitted heavy holes, 
light holes, and spin-orbit splitted bands. The summa- 
tion runs over all eigenstates. 

One important quantity that determines the nonequi- 
librium spin density fEg. pT]) ) is the scattering rate of 
hole carriers by Mn ions. We use Fermi's golden rule to 
calculate the ratci^ 
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where (pk^k' is the angle between two wave vectors k and 
k' . The matrix element M^'^, between two eigenstates 
(fe,n) and (fe',n') \fM 
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Here e is the d ielectr ic constant of the host semiconduc- 
tors and p — -y/ e^g/ e is the Thomas- Fermi screening wave 
vector,— where g is the density of states at Fermi level. 
Equation (ITi|) suggests there are two scattering chan- 
nels contributing to the scattering rate, i.e., the spin- 
dependent part that is due to the exchange coupling, 
thus proportional to Jp^, and a spin- independent part 
that arises from the Coulomb interaction. 

We calculate the field-like spin torque by using 
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where Jcx = J-pdNMnSa- Technically, by diagonalizing 
the total Hamiltonian Eq. ([T0| , we obtain the eigenvalues 
and eigenvectors which are then used to calculate the 
average values as well as the spin density. A formulation 
that is similar to the above has been developed to explain 
anisotropic magnetoresistance22iM and analyze magneto- 
optical properties in the same systemj2^ In this article, 
Eqs. (fTT|) - p5l) allow us to study the spin-orbit torques for 
a wide range of material parameters. 

At this stage, we need to emphasize that the inter- 
band transitions, arising from distortion in distribution 
function induced by an applied field, are neglected in our 
calculation. We believe such an omission is justified for 
two physical reasons. First of all, in the first-order Born 
approximation considered here, the intraband scattering 
is dominant.— Second, in a relatively clean system where 
the energy level broadening induced by impurity scatter- 
ing is much weaker than the Fermi energy, the interband 
contribution is considered insignificant 



III. FERMI SURFACE 

In this section, we concentrate on the behavior of the 
spin torque as a function of the Luttinger parameters (7) 
that determine the band structures as well as the Fermi 
sphere anisotropy originate Hamiltonian Eq.([T|). In the 
following, we first focus on the so-called spherical ap- 
proximation to reveal the role of the difference between 
majority and minority hole densities. Then we relax the 
spherical approximation to the most general case char- 
acterized by the Fermi surfaces that deviate significantly 
from a sphere. 



A. Spherical Fermi surface 

The spherical approximation in this section refers to 
the particular case 72 =73. We plot in FiglTJa), for 
different values of 72, the torque as a function of the 
magnetization orientation. The initial magnetization di- 
rection is along the z direction and the applied electric 
field is applied along the x direction. In such a setup, 
the effective Dresselhaus field {B-q) is (anti)aligned to 
the transport direction of the current. We expect that, 
in the leading order in /3, the noncquilibrium spin den- 
sity is mostly along the x direction too. Therefore, the 
y-component of the spin torque Ty is dominant since we 
consider in this article the magnetization direction ro- 
tates in the xz-plane. 

In FigUKa), the torque {Ty) is at its maximum when 
the azimuthal angle 9 — (i.e., Bd is perpendicular to 
magnetization direction) and decreasing monotonically 
towards its minimum at 6 = 7r/2 (i.e., Bd is aligned 
with the magnetization). The spin torque fits well by 
the curve T (x cos 9 with a maximal value that increases 
by increasing 72, when the other parameters are fixed. 

We explain this by plotting, in Fig. [Hb)-(d), the Fermi 
surfaces of the majority and minority hole bands in the 
kx ~ ky plane. As 72 becomes larger, the Fermi radii 
of the spin-splitted heavy hole band increases, so is the 
noncquilibrium spin density that is proportional to the 
population difference between the majority and minority 
bands. The result is an enhanced spin density, thus a 
larger torque. 



B. Non-spherical Fermi surface 

In the most general case 72 7^ 73, we have a 
non-spherical Fermi surface with the effect of band 
warping. 3^ An interesting correlation between the Fermi 
surface anisotropy and spin torque emerges when we plot 
the torque as a function of the magnetization angles. 
We shall be reminded that 73 is associated with the 
anisotropy energy at the gamma point The larger the 
difference between 73 and 72 , the more the Fermi surface 
deviates from a sphere. 
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FIG. 1. (Color online). The spherical Fermi surface approximation, i.e., 72 = 73 7^ 0. Panel (a) The y-component of the 
spin torque as a function of the magnetization direction for various 72(73) values. 9 is the angle between the z axis and the 
magnetization direction. Angle 4> is between the x axis and the projection of magnetization in the xy-plane. In this figure, 
angle 4> = Q. Fermi surface intersection at the fc^ = plane for (b)72 = 73 = 1.0, (c)72 = 73 = 2.0, (d)72 = 73 = 3.0. The red 
(empty) square, black dots, orange triangle and blue dots correspond to majority heavy hole, minority heavy hole, majority 
light hole and minority light hole band, respectively. The others parameters: 71 = 6.98, x = 0.06, Jpd ~ 55 meV nm^ and 
p = 0.2 nm"^. 



For a given 72 ~ 2.0 and various values of 73, we plot 
the spin torque as a function of magnetization angle in 
Fig. IHa). In this figure, as the Fermi surface deviates 
from a sphere (i.e., curves corresponding to 73 — 1.0 
and 73 — 2.93), the torque shifts avifay from a simple 
cos 9 function. In a comparison to the spherical case, 
the maximal value of the torque at = is lovifcred as 
73 7^ 72. As Eg. dTlT) indicates, in the linear response 
treatment formulated in this article, the magnitude of 
the spin torque is determined by the transport scatter- 
ing time and the expectation values of spin and veloc- 
ity operators of holes. Qualitatively, as the Fermi sur- 
face deviates from a sphere, the expectation value {sx) of 
the heavy hole band, contributing the most to the spin 
torque, is lowered at = 0. 

We also show in Fig. a more interesting behavior 
of the spin torque: as the Fermi surface warps, the an- 
gular dependence of the spin torque shifts away from a 
simple cos 9 dependence and develops, in an addition to 
the cosO envelop function, an oscillation with a period 
that is shorter than tt. The period of these additional 
oscillation increases as the Fermi surface becomes more 



anisotropic in fc-space, see Fig. [Hb) and (c). To fur- 
ther reveal the effect of band warping on spin torque, we 
plot Ty/ cos 9 as a function of the magnetization angle 
in Figl^fd). When 73 = 2.0 (spherical Fermi sphere), 
Ty/ cos 9 is a constant, for T oc cos9. When 73 = 2.93 
or 1.0, the transport scattering time of the hole carriers 
starts to develop an oscillating behavior in 6*^24 which 
eventually contributes to the additional oscillation in the 
spin torque. 

The correlation between the angular dependence of 
spin torque and the Fermi surface anisotropy found here 
is similar to the so-called crystalline anisotropic magne- 
toresistance that emerges from band warping.— In fact, 
such a correlation partially allows a possible mapping be- 
tween the symmetry properties of Fermi surface and the 
spin torque that is dynamically measurable. 

Material-wise, (In,Mn)As is popular in experiments 
and device fabricatioui^"— Although (In,Mn)As is, in 
terms of exchange coupling and magnetic properties in 
general, rather similar to (Ga,Mn)As, the difference in 
band structures, lattice constants, and Fermi energies 
between these two materials gives rise to different den- 
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FIG. 2. (Color online) (a)The j/-component of the spin torque as a function of magnetization direction. Fermi surface intersection 
in the — plane for (b)73 — 1.0 and (0)73 = 2.93. Panel (d) depicts Ty/ cos 9 as a function of magnetization direction. 
The red (empty) square, black dots, orange triangle and blue dots correspond to majority heavy hole, minority heavy hole, 
majority light hole and minority light hole, respectively. The others parameters are (71,72) = (6.98,2.0), Jpd ~ 55 meV nm^ 
and p = 0.2 nm^'^. 



sity of states, strains, and transport scattering rates. We 
plot in Fig 12] the spin torque (Ty) as functions of 9 for 
both (Ga,Mn)As and (In,Mn)As. The following param- 
eters are used in the calculation: the spin-orbit coupling 
parameter /? is found to vary, from 1 — 5 meV nm;'^* 
the exchange coupling constant Jpd — 55 meV nm^ for 
(Ga,Mn)As^ and 39 meV nm^ for (In,Mn)As4a For both 
materials, the spin torque decrease monotonically as the 
angle 9 increases from to tt/2. Throughout the en- 
tire angle range [0,7r], the amplitude of the torque in 
(In,Mn)As is larger than that in (Ga,Mn)As. We at- 
tribute this to two effects: (1) the spin-orbit coupling 
constant /3 in (In,Mn)As is about twice as larger than 
that in (Ga,Mn)As; (2) for the same hole concentration, 
the Fermi energy of (In,Mn)As is higher than that of 
(Ga,Mn)As. 



IV. EXCHANGE INTERACTION 

In our model, the spin torque depends not only on 
the properties of Fermi surface but also on the exchange 
coupling. From this aspect, our results are in agreement 
with the findings in Refl24l. 



Within the linear response formalism, i.e., Ea. ([TT|) . the 
p — d type exchange coupling makes two contributions to 
the spin torque. First, the exchange coupling is one of 
the factors that determine the energy eigenvalue of the 
total Hamiltonian, separating the majority and minority 
spin bands and contributing to the expectation value of 
the spin operator, and thus to the spin torque. The larger 
Jpd becomes, the larger is the energy splitting between 
majority and minority holes. Second, the exchange cou- 
pling contributes to the spin-dependent part of the trans- 
port scattering rate (time). A large Jpd leads to a short 
transport time. Therefore, we expect the competition 
between such two contributions eventually gives rise to 
a spin torque that depends nonlinearly on the exchange 
coupling Jpd. 

In Fig. Ufa), Ty component of the spin torque is plotted 
as a function of the exchange coupling Jpd, for different 
values of /3. In the weak exchange coupling regime, the 
electric generation of nonequilibrium spin density domi- 
nates, then the leading role of exchange coupling is de- 
fined by its contribution to the transport scattering rate. 
We provide a simple qualitative explanation on such a pe- 
culiar Jpd dependence. Using a Born approximation, the 
scattering rate due to the p—d interaction is proportional 
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FIG. 3. (Color online) Torque Ty as a function of the magneti- 
zation direction for (Ga,Mn)As (black square) and (In,Mn)As 
(red dots). For (Ga,Mn)As, (71,72,73) = (6.98,2.0,2.93); 
for (In,Mn)As, (71,72,73) = (20.0,8.5,9.2). The strength 
of the spin-orbit coupling constant is: for (Ga,Mn)As, /3 = 
1.6 meV nm; for (In,Mn)As, P — 3.3 meV nm. 



to 1/tj = &</pd: where parameter b is Jpd- independent. 
When the nonmagnetic scattering rate 1 /tq is taken into 
account, i.e., the Coulomb interaction part in Ea. (ll4p . 
the total scattering time in Eg. lfTTj) can be estimated as 
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which contributes to the torque by T cx Jpd/(?ir). This 
explains the transition behavior, i.e., increases linearly 
then decreases, in the not-so-large Jpd regime in FigH) 

As the exchange coupling further increases, Ea. ([T5)) is 
dominated by the spin-dependent scattering, therefore 
the scattering time 1/hT (x 1/Jpd- Meanwhile, the en- 
ergy splitting due to exchange coupling becomes signifi- 
cant, thus (s) cx Jpd. In total, the spin torque is insensi- 
tive to Jpd, explaining the flat curve in the large exchange 
coupling regime. 

In Fig. mb) , we plot the influence of the exchange cou- 
pling on the spin torque for two materials. In (In,Mn)As, 
mainly due to a larger Fermi energy in a comparison to 
(Ga,Mn)As, the peak of the spin torque shifts towards 
a larger Jpd. The amplitude of torque in (In,Mn)As is 
clearly larger, arising from a more robust spin-orbit cou- 
pling. 



V. HOLE CARRIER CONCENTRATION 

The possibility to engineer the electronic properties by 
doping is one of the defining features that make DMS 
promising for applications. Here we focus on the doping 
effect which allows the spin torque to vary as a function 
of hole carrier concentration. 

In Fig. [Sfa) , the torque is plotted as a function of the 
hole concentration for different /3 parameters. With the 



FIG. 4. (Color online) The Ty component of the spin torque 
as a function of exchange coupling Jpd. (a) Ty versus Jpd 
at various values of /3, for (Ga,Mn)As. (b) Ty versus Jpd, 
for both (Ga,Mn)As and (In,Mn)As. The magnetization is 
directed along the z-axis {9 = Q). The other parameters are 
the same as those in FigO 



increase of the hole concentration, the torques increase 
due to the enhanced Fermi energy. At a small /3 (weak 
spin-orbit coupling) , the torque as a function of the hole 
concentration {p) follows roughly the p^/^ curve as shown 
in the inset in Fig. [Slja). The spherical Fermi sphere 
approximation and a simple parabolic dispersion relation 
allow an analytical expression for the spin torque, i.e., in 
the leading order in /3 and Jpd, 



m* /3 Jpd 



(17) 



where to* is the effective mass. The Fermi energy Ep 
and the Drude conductivity are given by 



2m 



2/3 



(18) 



where r is the transport time. The last two relation im- 
mediately gives rise to T oc p^/'^. In the six-band model, 
the Fermi surface deviates from a sphere and as the value 
of /3 increases, the spin-orbit coupling starts to modify 
the density of state. Both effects render the torque versus 
hole concentration curve away from the p^^^ dependence. 
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In Fig. [UJb), we depict the torque as a function of hole 
concentration for mentioned materials. In (In,Mn)As the 
amplitude of the torque is larger than that of (Ga,Mn)As, 
given the same hole concentration. We attribute this 
mainly to a larger Fermi energy in (In,Mn)As. 
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FIG. 5. (Color online) The y-component of the spin torque 
as a function of hole concentration, (a) The y-component 
of the spin torque versus hole concentration at different 
p. (b) spin torque versus hole concentration in (Ga,Mn)As 
and (In,Mn)As. For (Ga,Mn)As, Jpd = 55 meV nm^; for 
(In,Mn)As, Jpd ~ 39 meV nm^. The other parameters are 
the same as in FigU) 



VI. SUMMARY AND CONCLUSION 



Mediated by itinerant charge carriers, the interplay be- 
tween a spin-orbit coupling and an exchange splitting 
gives rise to a spin torque. Using a Kohn-Luttinger 
Hamiltonian supplemented by a linear Dresselhaus spin- 
orbit coupling and an exchange interaction, in the frame- 
work of linear response, we have studied theoretically the 
spin torque in two typical III-V diluted magnetic semi- 
conductors (Ga,Mn)As and (In,Mn)As. 

Beyond the spherical Fermi sphere approximation that 
was commonly studied in literature, we found that the 
spin torque as a function of the magnetization direction 
has an intriguing correlation with the anisotropy of the 
Fermi surface. Our calculation also suggests that the 
spin torque develops a nonlinear dependence on the ex- 
change coupling, which is the result that the exchange 
coupling not only determines the energy eigenstate but 
also contributes to the transport scattering rate of the 
hole carriers. We have also investigated the spin torque 
as a function of hole concentration. From aspect of mate- 
rial selection, given other parameters the same, the spin 
torque in (In,Mn)As is usually larger than its counterpart 
in (Ga,Mn)As. 
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